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ABSTRACT 

Monte Carlo calculations of the Ka line fluoresced by coronal X-rays and 
emitted just above the temperature minimum region of the solar atmosphere have 
been employed to investigate the use of this feature as an abundance diagnos- 
tic. While quite weak, we estimate line equivalent widths in the range 0.02-0.2 A, 
depending on the X-ray plasma temperature. The line remains essentially uncon- 
taminated by blends for coronal temperatures T < 3 x 10 6 K and should be quite 
observable, with a flux > 2 ph s _1 arcmin -2 . Model calculations for solar chemical 
mixtures with an O abundance adjusted up and down by a factor of 2 indicate 
35-60% changes in O Ka line equivalent width, providing a potentially useful O 
abundance diagnostic. Sensitivity of equivalent width to differences between re- 
cently recommended chemical compositions with "high" and "low" complements 
of the CNO trio important for interpreting helioseismological observations is less 
accute, amounting to 20-26% at coronal temperatures T < 2 x 10 6 K. While still 
feasible for discriminating between these two mixtures, uncertainties in measured 
line equivalent widths and in the models used for interpretation would need to 
be significantly less than 20%. Provided a sensitive X-ray spectrometer with 
resolving power > 1000 and suitably well-behaved instrumental profile can be 
built, X-ray fluorescence presents a viable means for resolving the solar "oxygen 
crisis" . 



Subject headings: Sun: abundances — Sun: activity — Sun: corona — Sun: 
X-rays — X-rays: stars — radiation mechanisms: nonthermal 
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Introduction 



At the turn of the millennium, precise agreement between the observed oscillation spec- 
trum of the Sun and pre dictions of theoretical sol ar models, built using the best abundance 
assessments o f the day (IGreyesse fc Sauvallll998l ). represented a triumph of modern astro- 



physics (e.g. iBahcall et al.l l2005al ). This illustrious accord has recently been clouded by 



sophisticated re-assessments of the solar complement of light elements based on 3-D non- 
LTE hydrodynamic photospheric modelling: such models demand less of the elements C, 
N, O and Ne that are important for the op acity of the solar interior by 25-35 % compared 
to earlier assessments ( Asplund et al.ll2005l ). Solar models employing this chemical compo- 
sition lead to predictions of the depth of the convection zone, helium abun dance, density 
and sound speed in serious disagreement with helios eismo logy measurements (IBasu &: Antia 
2004 ; iTurck-Chieze et al.l 12004] ; IBahcall et al.ll2005al ). This state of affairs is no w commonly 



referr ed to as the "solar model problem", or the "solar oxygen crisis" 
2006h . 



e.g. 



Ayres et al. 



To solve this problem, lAntia fc Basul (120051 ) and IBahcall et al.l (j2005bl ) suggested the 
uncertai n solar Ne abundance might be raised to compensate for lower C, N and O abun- 
dances. iDrake fc Testal (120051 ) found empirical support for this in Chandra high resolution 
X-ray spectra of mostly magnetically active stars for which the Ne/O abundance ratio ap- 
pears consistently higher by a factor of ~ 2 or more than the currently recommended solar 
value of Ne/O=0.15 by number. Th ey argued that existing solar coronal Ne/O measure- 
ment s supporting the latter (see, e.g., brake fc Testal 120051 : ISchmelz et aDbood : lLandi et al.l 
20071 ) might be biased by chemical fractionation. However, subsequent studies of the he- 



lioseismo logy data using different theoretical models disag r ee as to the vi ability of such a 
solution (IDelahaye fc Pinsonneaultll2006l ; IZaatri et al.l 120071 : iLin et al.l 120071 ). 



An alternative possi bility is that the su ggested do wnward C, N and O abun dance re- 
visions are unwarranted. I Ayres et al.l (120061 ) favour the iGrevesse fc Sauvall (119981 ) value of 
0/H= 8.83 based on multi-component atmospheric modelling and analysis of CO lines. How- 
ever, a very recent 3D e mpirical study of spatially- r esolve d spectropolarimetric observation 
of Fe I and O I lines by ISocas-Navarro fc Norton! (120071 ) finds 0/H= 8.63 ± 0.1, in good 
agreement with Asplund et al. ( 20051 ) . 



The difficulty in determining the solar O abundance lies in the complexity of the non- 
local thermodynamic equilibrium optically-thick line formation problem in the inhomoge- 
neous and turbule n t pho tosphere. In this paper, we study a possible alternative approach. 
Drake &: Ercoland (120071 ) recently investigated the utility of the weak photospheric Ka flu- 
orescence line of neutral Ne excited by coronal X-ray photoionisation as an independent 
abundance diagnostic. Here, we turn our attention to the O Ka line at 23.37 A (530.6 keV). 
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In contrast to the complex photospheric absorption line formation problem, the formation 
of the X-ray fluorescent O Ka line is relatively simple, depending only on the geometry of 
the overlying coronal s ource of irradiation and the heliocentric angle of the photospheric 
region in question (e.g. lBailll979l ; iDrake et al.l 120071 ; iDrake fc Ercolanoll2007l ). We present 
Monte Carlo calculations of the fluorescent O Ka line in §2, and address its potential as an 
abundance diagnostic in §3 and §4. 



2. Monte Carlo Calculations of Oxygen Fluorescence 



We adopt the same computational metho ds described i n ear lier work on Fe and Ne 
fluorescence lines, and t he reader is referred to IDrake et al.l (120071 . hereafter Paper I) and 
Drake fe Ercolanol (120071 . Paper II) for details. A brief summary is provided below. 



X-ray "characteristic" (fluorescence) Ka lines correspond to the 2p — Is decay of the 
excited state resulting from ejection of an inner-shell Is electron in the neutral (or near- 
neutral) atom by either electron impact or photoionisation. In the c ase of solar and stellar 



photospheric fluore scent lines, the former is expected to dominate (jBaskol Il978l ; iBail 11979 



Parmar et al.l 1 1 9841 ) . For a given coronal X-ray source spectrum, F(X), the observed flux of 
Ka photons from the photosphere then depends on the photospheric abundance A of the 
fluorescing species relative to that of other elements of significance for the photoabsorption 
opacity in the vicinity of the Is ionisation edge; the height h of the emitting source; and the 
heliocentric angle 9 between the emitting source and the observer (IBail 119791 ). 



Fluorescent lines are formed in the region of the atmosphere corresponding to optical 
depth of approximately unity for the primary K-shell ionising photons. In the case of O, 
for no rmal incidence the solar atmospheric Model C (VALC) of IVernazza. Avrett. fe Loeser 
(Il98ll ) reaches the K-shell photoabsorption r = 1 depth slightly above the temperature 
minimum, at a gas temperature of about 5600 K and about 800 km above the point where 
the continuum optical depth at 5000 A, r 50 oo, is unity. 

The expected intensity of the emergent O Ka line was compu ted using a modified 



versio n of the 3D Monte Carlo radiative transfer code MOCASSIN ( lErcolano et al.l 12003 



20051 . see also Paper I for a specific description of fluorescence calculations). This code follows 
energy packets representing incident X-ray photons as they interact with the photospheric 
gas through photoelectric absorption or scattering events until the packets escape. O Ka 
packets are produced following K-shell X-ray absorption with a probability dictated by the 
Ka fluorescent yield for oxygen: these packets are followed in the same fashion. 



We adopted a value of 8.3 x 10 3 for the fluorescence yield of neutral oxygen (IKrause 
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19791 ). Coronal spectra used for the incid ent X-ray radiatio n field were computed using 



emissivities from the CHIANTI database (jLandi et al. 



Mazzotta et al.l (119981 ) , as implemented in the PINTofALEU IDL suite of programs (jKashyap fc Drake! 
2000h. 



20061 ) and the ion pop ulations of 



For a given coronal X-ray spectrum, the maximum intensity of a fluorescent line is 
achieved for a heliocentric angle 9 = and height h = 0. As in the case of the earlier study 
of Ne, since we are primarily interested here in the observability of the fluorescent line, we 
adopt 9 = and h = as baseline parameters, these values yielding the maximum possible 
line strength. 

The wavelength of the O Ka doublet has not been determined with high precision and 
we adopted the value 23.37±0.02 A, corresponding to the mean of values found for th e energy 



l evel o f the excited \ls ] 2p 44 P of 5 30.6 ± 0.3 eV from laboratory measuremen ts by iKrause 



( 1994), ICaldwell et al.l (119941 ) and IStolte et al.l (119971 . see also the review by iGarcia et al 



20051 ). This wavelength is in good agreement with the position of the 2p 34 S —> [ls]2p 44 P 
resonance seen in abs o rption in the interstellar medium toward bright X-ray continuum 
sources by lJuett et al.l (120041 ) . 



The O Ka flux was computed for isothermal irradiating coronal spectra with plasma 
te mperatures in the range 1 6 -10 7 K and coronal and photospheric chemical compositions 
of iGrevesse fc Sauvall (119981 . hereafter GS; with 0/H=8.83 on the usual log+12 scale) and 
Asplund et al.l (120051 . 0/H=8.66). As noted in §1, these two chemical mixtures differ not 
only in O abundance but also, primarily, in C, N and Ne abundances: these differences 
turn out to be important for the O fluorescent line and are discussed in more detail below. 
To investigate changes in the O abundance only, we also performed calculations using a 
photospheric composition corresponding to the GS abundance mixture with O elevated and 
decreased by a factor of 2 (0/H=9. 13,8.53, or [O/H]=±0.3). 



3. Strength of the O Ka Line 

3.1. Monte Carlo Results 

The photospheric O Ka fluorescent line is shown in comparison with the direct coronal 
spectrum for the case of GS abundances for both coronal spectrum and photosphere in 



1 The Package for INTeractive Analysis of Line Emission is freely available from 



http://hea-www.harvard.edu/PINTofALE/ 
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Figured] for the range of coronal plasma temperatures 6.0 < logT < 7.2. The lower panel 
illustrates the same spectra seen at a resolving power of A/AA = 1000, where AA is assumed 
to be the full-width at half-maximum of a Gaussian instrument response function. Unlike 
the Ne Ka line, the O Ka transition forms in a relatively uncrowded spectral region and 
the CHIANTI database predicts no devastatingly large blends in the immediate vicinity of 
23.37 A. For temperatures approaching 10 7 K, a line of Fe XXIII at 23.36 A appears, but 
for GS abundances is always dominated by the O Ka line by a factor of 10. Lying redward 
of our adopted line centre is Ca XV A23.39, which peaks at slightly cooler temperatures 
(logT ~ 6.6) at an intensity of about 30 % of that of O Ka. 

The full list of lines in the CHIANTI database within 5cr of the O Ka wavelength 
range (±0.05 A from 23.37 ±0.02 A) with intensities > 10~ 4 times that of the brightest line 
(Fe XXIII A23.363) ar e listed in Table HI The NIST Atomic Spectra Database (version 3.1.2; 



Ralchenko et al.l 120071 ) also lists two other transitions within the 5a range from Ti XIII and 



Sc XIV. Neither of these are expected to be of any significance since their solar abundances 
according to the GS assessment are 300 and 21000 lower than that of Fe, respectively. 

The equivalent width (EW) of the O Ka line relative to coronal fluorescing spectra with 
GS abundances is illustrated in Figure [2] for a range of isothermal plasma temperatures and 
different photospheric O abundances. The Monte Carlo sampling error on the computed 
EWs is estimated to be no larger than 7%. One striking feature of the trend of EW with 
T for all compositions is the sharp decline with rising temperature for logT < 6.6: this 
contrasts with both Ne Ka and Fe Ka, whose EWs exhibit a steady, monotonic rise with 
increasing T (Papers I & II), corresponding to a commensurate increase in the number 
of ionising photons. The different behaviour of O Ka arises because, for coronal plasmas 
with logT < 6.6 and solar composition, K-shell photoionisation of oxygen is due mostly 
to line rather than continuum radiation. Toward hotter temperatures, the important line- 
emitting photoionising species, including O VII and O VIII, become ionised, and continuum 
contributions begin to dominate. 

Regarding sensitivity to O abundance, Figure [2] illustrates that the EW changes by a 
smaller factor than expected base d on a proportional relation with the photospheric the O 
abundance. As discussed by iBail ( 119791 ) and in Paper I, this is a result of the the O K- 
shell photoionisation cross-section being a significant component of the total opacity near 
threshold: for very large O abundances where O begins to dominate the opacity, the EW 
will tend to a constant value dictated simply by equipartition of ionising photons between O 
K- and L-shells. Nevertheless, for "normal" ranges of O abundance the O Ka line variations 
are significant: we find typical changes of 35-60% over the full range of coronal temperatures 
for O abundance variations of a factor of 2. 
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When the other abundant light elements C and N are allowed to scale with the 
abundance the diagnostic fares slightly less well but does pro vide di scrimination be t ween the 
photospheric compositions of GS and lAsplund et al.l (120051 ). The lAsplund et al.l (120051 ) O 
abundance is 48% lower than that of GS, and the O Ka EWs differ by 20-25%. For a given 
exciting X-ray spectrum, changes in fluorescent line strength with different photospheric 
parameters depend on changes in the fraction of ionising photons that ionise the parent shell 
of the fluorescent line in question. For a solar composition, two other sourc es of opacity in 



the vi cinity of the O edge are C and N. The chemical compositions of GS and lAsplund et al. 



( 120051 ) differ by < 10% in O/C and O/N ratios, and the lockstep changes in these elements 
that dilutes the effect of the different O abundances relative to H on the O Ka EW. Despite 
this slightly lower sensitivity of the fluorescent line to the global chemical composition, 
accurate measurements of the O Ka EWs still po tentially provide a n ew and relatively 
direct means of assessing the veracity of the GS and lAsplund et al.l (120051 ) mixtures. 



The line EW is also sensitive to the abundances adopted for the exciting coronal spec- 
trum. This is due to the large contribution of the O VII He-like complex to the source of 
ionising photons for coronal temperatures logT < 6.5. To a lesser extent, Fe L-shell and 
Ne, Mg and Si H-like and He-like lines also make a contribution for temperatures up to 
logT ~ 7.0. We have examined the sensitivity to coronal abundan ces by compar i son o f 
O Ka line EWs computed for coronal spectra generated using GS and lAsplund et al.l (120051 ) 
compositions, as illustrated in Figure [2j The former are higher than the latter by an amount 
that decreases from 25-30% for temperatures logT < 6.3 where the O VII lines dominate, 
to ~ 10% at logT ~ 6.8-7.0. Differenc es at higher tempera tures are largely due to the lower 
Ne, Mg, Si and Fe abundances in the (lAsplund et al.l 120051 ) composition (by 74, 12, 10 and 
55% respectively). 

In addition to uncertainties in the solar O content, coronal abundance variations are 
also expected as a result of chemical fractionation, in which the abundances of elements with 
low first i onisation poten tials are seen to differ from photospheric values by factors of up to 



4 (e.g. lFeldmanlll992l ). In this context, we emphasise that in a practical application of 



the fluorescence technique the photospheric abundance would be deduced by comparing the 
model fluorescent EWs computed for different photospheric abundances and the simultane- 
ously observed coronal spectrum, and knowledge of the actual coronal composition is not 
required. 
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3.2. Observability of O Ka Fluorescence and its Abundance Diagnostic Utility 

As discussed in Paper II in the context of Ne fluorescence, existing solar spectra generally 
have insufficient sensitivity to detect these relatively weak lines. While the Ka line is 
expected to be an order of magnitude stronger than that of Ne Ka, there are essentially 
no high quality solar spectra that reach beyond 23 A, including the most extensive sets of 
soft X-ray spectra obt ained for the Sun b y the satellite-borne Solar Maximum Mission Soft 
X-ray Polychromator ( Acton et al.l ll98ol ) and SOLEX instruments ( McKenzie et al. 1980 ) 



that had coverage up to 22.43 and 23.0 A, respectively. 

Existing analyses of higher quality spec tra whose ranges encompassed the O K a line 



obtained with rocket-borne instruments (e.g. iMcKenzie et al.lll978l ; lActon et al.lll985l ) make 



no mention of a feature at the expected wavelength of 23.37 A. The w eakest intens i ty rep orted 



from the spectrum of a flare in the 10-100 A range obtained by lActon et al.l ( 19851 ) was 
10 ph arcsec 2 , which can be compared with the measured O VII resonance line intensity of 
508 ph arcsec 2 . At a typical solar coronal temperature of logT = 6.3, the flux of the O VII 
A21.6 line is a factor ~ 10 3 stronger than our prediction for the O Ka line and the non- 
detection of the O Ka line is not a surprise. While the intensity ratio drops with increasing 
isothermal temperature to a factor of 10 at logT = 6.8, where the O VII ion fraction is 
comparatively small, it would take the presence of relatively little cooler material within the 
instrument field of view to strengthen the observed O VII A21.6 line significantly. Accurate 
measurement of the O Ka line requires greater sensitivity than afforded by existing solar 
spectra. 

As shown in §3.11 the O Ka line remains essentially blend-free and isolated for plasma 
temperatures encountered in the non-flaring solar corona. The line EW peaks at coronal 
temperatures of logT ~ 6.2 — conditions typical of those found in the coolest active regions 
and in the quiet Sun. In the case of Ne Ka, we found fluorescent line fluxes for quiet 
Sun conditions of ~ 0.2 ph s _1 arcmin -2 . The O Ka line is expected to be an order of 
magnitude stronger than this. Moreover, since the Fe XXIII and Ca XV blends remain at 
levels significantly below 10% of O Ka for temperatures up to logT ~ 6.5, the line is also 
expected to be accessible in much brighter active regions where fluxes are higher. 

As in the case of Ne Ka, the large intensity contrast between O Ka and neighbouring 
strong lines — in this case the O VII He-like complex in the range 21.6-22.1 A, and the N VII 
A24.78 resonance line — imposes tight instrumental profile requirements in order to avoid 
contamination of the O Ka region from surrounding line wings. In the absence of suitable 
simultaneous spectroscopic observations of the solar regions being studied for fluorescence, 
such an instrument would also require fairly broad wavelength coverage so as to measure 
accurately the fluorescing X-ray spectrum from the O K threshold down to wavelengths at 
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which the contribution to O K ionisation is no longer significant. The short wavelength 
cut-off requirement will depend on the coronal temperatures being studied. 

We emphasise that in a practical application of the fluorescence technique for determin- 
ing photospheric abundances, knowledge of the actual coronal composition is not required 
since the observed EW is compared to a measured coronal spectrum that itself must be used 
in model predictions. 



4. Conclusions 

We have investigated the observability of the O Ka line formed in the solar upper 
photosphere/lower chromosphere by fluorescence from overlying coronal X-radiation. While 
quite faint compared to the prominent lines of abundant coronal ions, O Ka should be 
essentially unblended at plasma temperatures up to logT ~ 6.5. 

The line EW peaks at temperatures logT ~ 6.2 with a value of ~ 0.1 A, and we 
estimate an O Ka flux of ~ 2 ph s _1 arcmin -2 for the quiet Sun at heliocentric angles close 
to 0° — an order of magnitude larger than predicted in Paper II for Ne Ka fluorescence. Both 
lines should be quite observable with instruments adequately tailored toward high resolution 
A/AA > 1000 and sufficiently well-behaved instrumental profile wings. 

The O Ka line provides a potentially accurate absolute O abundance diagnostic, pro- 
vided abundances of other prominent light elements are also constrained. It should also 
be possible to distinguish between the currently competing photospheric compositions of GS 



and lAsplund et al.l (120051 ). provided line EWs can be determined to an accuracy significantly 
better than 20%. We conclude that X-ray fluorescence of O (and Ne) represents a viable 
means for resolving the current solar abundance-helioseismology conundrum. 
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Table 1: Lines in the CHIANTI 5 database lying within ±0.07 A (5er + a rest wavelength 
uncertainty 0.02 A) of the position of the 23.37 A O Ka transition for a Gaussian resolving 
power of A/AA = 1000. Only lines with predicted intensities > 1 x 10~ 4 that of the brightest 



line are listed. 


Ion 


A (A) 


Rel. Int. 


log T max 


Transition 


Fe XXIII 


23.303 


0.54 




7.2 


2s3p 3 P 2 - 2shd 3 D 3 


Ca XV 


23.305 


0.07 




6.6 


2s 2 2p 2 1 D 2 - 2s 2 2p 3 d 3 D 1 


Ca XVI 


23.320 


0.03 




6.7 


2s2p 2 2 P 3/2 - 2s2p( 3 P)3d 4 D 5/2 


Fe XXIII 


23.321 


0.10 




7.2 


2s3p 3 P 2 - 2shd 3 D 2 


Ca XVI 


23.341 


1 x 10~ 


-3 


6.7 


2s2p 2 2 P 3/2 - 2s2p( 3 P)3d 4 D 3/2 


Ar XVI 


23.343 


0.01 




6.7 


ls 2 2s 2 S 1/2 - ls 2 3d 2 D 5/2 


Ar XVI 


23.362 


0.01 




6.7 


ls 2 2s 2 S 1/2 - ls 2 3d 2 D 3/2 


Fe XXIII 


23.363 


1.00 




7.2 


2s3p 3 Pi - 2s5s 1 S 


Fe XXIII 


23.381 


0.04 




7.2 


2s3p 1 P l - 2s5d 3 Dr 


Fe XXIII 


23.388 


2 x 10" 


-4 


7.2 


2s3p 1 P 1 - 2shp 3 Pi 


Fe XXIII 


23.388 


3 x 10" 


-4 


7.2 


2s3p l Px - 2shp l P x 


Ca XV 


23.388 


0.88 




6.6 


2s 2 2p 2 1 S - 2s 2 2p 3 d 1 P 1 


Fe XXIII 


23.406 


8 x 10" 


-3 


7.2 


2s3p 3 P 2 - 2shd 3 D 1 


Fe XXIII 


23.412 


3 x 10" 


-4 


7.2 


2s3p 3 P 2 - 2shp 3 Pi 


Ca XV 


23.415 


0.45 




6.6 


2s 2 2p 2 1 D 2 - 2s 2 2p 3 d 3 F 3 


Ni XXVIII 


23.439 


5 x 10- 


-4 


8.0 


3p 2 P 1/2 4d 2 D 3/2 
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23.25 23.30 23.35 23.40 23.45 




23.25 23.30 23.35 23.40 23.45 

Wavelength [A] 



Fig. 1. — Top: The strength of the O Ka shown in comparison to neighbouring and blending 
lines in the fluorescing coronal X-ray spectrum for coronal isothermal plasma temperatures 
in the range 10 6 -10 72 K (wavelengths shown off-set for clarity). The underlying grey shaded 
region represents the uncertainty in the wavelength of the O Ka line. Bottom: The same 
coronal spectra with (solid curves) and without (dashed) the addition of O Ka smoothed to 
a resolving power (FWHM) of A/AA = 1000. 
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Fig. 2. — The equivalent width of the O Ka line with respect to the ionising coronal X- 
ray spectrum for heliocentric angle 9 = as a function of isothermal plasma temperature, 
computed for different coronal and photospheric chemical compositions: GS corona + GS 
photosphere (solid black curve), GS corona + G S photosphere with t he O abundance dif- 
fering by factors of 2 (solid grey), GS coro na + lAsplund et al.l (j2005l ) photosphere (heavy 
dashed), GS coron a + GS photosphere w ith lAsplund et al.l (120051) O abu ndance (0/H=8.66; 
light dashed) ; and lAsplund et al.l (120051 ) corona + lAsplund et al.l (120051 ) photosphere (black 
dotted). 



